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We present two millisecond pulsar discoveries from the PALFA survey of the 
Galactic plane with the Arecibo telescope. PSR J1955+2527 is an isolated pulsar 
with a period of 4.87 ms, and PSR J1949+3106 has a period of 13.14 ms and is 
in a 1.9-day binary system with a massive companion. Their timing solutions, 
based on 4 years of timing measurements with the Arecibo, Green Bank, Nangay 
and Jodrell Bank telescopes, allow precise determination of spin and astromet- 
ric parameters, including precise determinations of their proper motions. For 
PSR J1949+3106, we can clearly detect the Shapiro delay. From this we mea- 
sure the pulsar mass to be 1.47l ; 3 i M , the companion mass to be 0.85lo!n M 
and the orbital inclination to be i = 79.9+}g degrees, where uncertainties corre- 
spond to ±1-0" confidence levels. With continued timing, we expect to also be 
able to detect the advance of periastron for the J1949+3106 system. This ef- 
fect, combined with the Shapiro delay, will eventually provide very precise mass 
measurements for this system and a test of general relativity. 



1. Introduction 

In this paper, we discuss timing results for two millisecond pulsars (MSPs) discovered 
by the PALFA Consort iurrl^l with the Arecibo telescope. The PALFA survey of the Galactic 
plane and follow-up observations of new discoveries are motivated by the wide applications 
of pulsar timing in exploring the composition, internal structure, and magnetospheric state 
of neutron stars. Millisecond pulsars in particular tend to be extremely stable rotators, 
which can be used to address a variety of problems in fundamental physics and astrophysics. 
One such outstanding problem is detecting gravitational waves and studying the properties 
of gravitational wave radiation from various types of sources and various epochs in the 
Universe's lifetime. Timing observations of MSPs in binary systems can be used to estimate 
the pulsar and companion masses via measuring post-Keplerian binary parameters. Neutron 
star mass measurements allow us to constrain the equation of state (EoS) of matter at 
densities larger than that of an atomic nucleus. 

Until recently, all precise measurements of neutron star masses fell within a narrow 
range around 1.4 M , the Chandrasekhar limit. However, recent, precise measurements of the 
masses of some MSPs sh owed that they can have significantly higher masses. PSR J 1903+0327 



( iChampion et al.l 120081 ) . the first MSP found in the PALFA survey, has a mass of 1.67 ± 



0.02 M (99.7 % confidence limit, iFreire et al.l 1201 if ) which is significantly above the Chan 



16 http: //www.naic . edu/alf a/pulsar 
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drasekhar limit. PSR J1614— 2230 was found to have a mass of 1.97±0.04 M (IDemorest et al. 



2010|). These high masses rule out many EoSs for matter at densities higher than that of 



the atomic nucleus. In particular, the mass measure ment for J1614— 2230 rules out or highly 



constrains most p r oposed hyperon or boson EoSs (iGlendenning fc Schaffner-Bielichl 11998 



Lackey et all 12009 . ISchulze et al.ll2006L lLattimer fc Prakashl 120071 ). 



In this paper, we describe two MSPs discovered b y the PALFA survey , which uses the 
Arecibo telescope and the seven-beam ALFA receiver (ICordes et al.ll2006l ). In addition to 
presenting full timing solutions for the two pulsars, we explore what intrinsic or extrinsic 
effects account for the overall TOA uncertainty in both cases. One of the two discoveries 
reported here is in a nearly edge-on binary system with a massive companion, allowing 
measurement of the Shapiro delay and consequently estimation of the pulsar and companion 
masses. 



2. Observations and Data Reduction 

PSR J1949+3106 and PSR J1955+2527 were disco vered in PALFA sur vey data taken 
in 2006 and processed by the Cornell search pipeline ( jPeneva et al.ll2009l ). The PALFA 
survey is ongoing since 2004 and uses the seven-beam Arecibo L-band Feed Array (ALFA) 
to search for pulsars in low Galactic latitudes (|6| < 5°) in the portions of the Galactic 
plane visible from Arecibo, 32° < I < 77° and 168° < I < 214°. PALFA has discov- 



ered 9 7 new pulsars, including 15 MSPs, as of Ju l y 2012 (ICordes et al. 
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Crawford et al.l 120121 ) . There are currently 23 pulsars in the ATNF cat- 



alog 19 ! ^Manchester et al.l l2005h with P < 15 ms and |6| < 5° that are not in globular 
clusters. Of those, eight have dispersion measur es (DMs) > 100 pc cm ~ 3 , and only one has 
DM > 200 pc cm -3 . That pulsar is J1903+0327 ( Champion et al.ll2008 ). an eccentric binary 
MSP discovered by the PALFA survey. The two discoveries we report, J1949+3106 and 
J1955+2527, have DM = 164 pc cm -3 and 209 pc cm -3 , respectively This makes the latter 
only the second Galactic plane MSP known with DM > 200 pc cm -3 . I Crawford et al.l (120121 ) 
present four new MSPs with DM > 100 pc cm -3 , also found by PALFA. One of these four 
pulsars has DM = 249 pc cm" 3 . Along with the high DMs of J1903+0327 and J1955+2527, 
this shows that the PALFA survey is probing deeper and finding MSPs whose detectability 
with earlier instruments may have been hampered by dispersion smearing. 



In 2006, the PALFA survey used the Wideband Arecibo Pulsar Processor (WAPP) 



'http : //www. atnf . csiro . au/research/pulsar/psrcat 
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backends (IDowd et al.ll2000l ) with 100 MHz ba ndwidth per beam, two summed polarization 
channels, 256 lags, and 64 /xs sampling time ( jCordes et al.l 120061 ) . Since 2009 the survey 
uses the newer PDEV/Mock spectrometei@, which allows data to be taken across the entire 
300 MHz bandwidth of the ALFA receiver for each beam, thus increasing sensitivity by a 
nominal factor of \/3. This observation setup uses two summed polarizations, 960 channels, 
and a samping time of 65.5 /xs. 

PSR J1949+3106 and PSR J1955+2527 were found in a conventional Fourier-based 
periodicity search of WAPP data and were not detected in a single pulse search. Timing 
observations of the pulsars began initially with the Green Bank Telescope. After a dedicated 
timing program was established, higher- sensitivity regular timing observations were done 
with the Arecibo, Nangay, and Jodrell Bank telescopes. 



2.1. Timing Observations 

Green Bank tim ing observations used the S-band receiver and the Pulsar Spigot backend 



flKaplan et all 120051 ) with a total bandwidth of 800 MHz centered on 1890 MHz, 1024 lags, 
and 81.92 /xs sampling time. Because of strong and persistent RFI, the lowest 200 MHz 
of the band is excluded, leaving 600 MHz of effective bandwidth centered on 1950 MHz. 
The Arecibo observations used the L-band Wide receiver and WAPP correlator backends. 
The four WAPPs covered adjacent 50 MHz bands with 512 lags each, centered at 1420, 
1470, 1520, and 1570 MHz. The sampling time for Arecibo observations was 64 /xs. GBT 
and Arecibo observations were made in search mode, with two summed polarizations and 
continuous data streams recorded for off-line processing. After data collection, the outputs 
of the spectral channels were folded at the predicted topocentric period of the pulsar, shifted 
to compensate for dispersive delays, and summed to produce a total intensity pulse profile. 

Observations with the Lovell telescope covered a 300-MHz band centered on 1500 MHz. 
A digital filterbank with a channel width of 0.5 MHz was used in a synchronous integration 
mode, in which the data acquisition is synchronous with the pulsar period. Baseband data 
were collected for the duration of a time bin (1/1024 of a pulse period). The power spectrum 
of these data is then calculated and added to the cumulative spectrum for that bin. After 
shifting to compensate for dispersive delays, the data were summed to produce a total 
intensity profile. This setup results in a sampling time of 12.83 /xs for J1949+3106 and 
4.76 /xs for J1955+2527. 



! http : / / www . naic . edu/phil/hardware/pdev/usersGuide . pdf 
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Nangay observations were performed with the Berkeley-Orleans-Nangay instrument (jCognard et al. 



20091 ) and used coherent dedispersion performed by Graphics Processing Units (GPUs) in 4 



MHz channels over a total bandwidth of 128 MHz centered at 1400 GHz. For each pulsar, 
all the observations at each frequency were integrated to produce a template used to derive 
topocentric t imes of arrival (TOAs) following a standard x 2 fit in the frequency domain 
rtTavlor! Il993 ). 



We used the TEMP02 software package ( iHobbs et al.l 120061 ) to perform least-squares 
fitting of various pulsar parameters by minimizing the square of the differences between 
expected and measured TOAs. We directed TEMP02 to fit for constant offsets between 
TOA sets from different observatories by bracketing each TOA set with JUMP statements. 
We scaled the TOA uncertainties (via EFAC statements, Table [T]) so that for each TOA 
set the ratio of \ 2 to the number of degrees of freedom is close to unity. TEMP02 applies 
clock corrections based on the location of each observatory and Earth rotation data to convert 
TOAs to terrestrial time Conversion from TT to coordinated b arycentric time (TCB) 

was done using the DE405 Solar System ephemeris (IStandishlll998r). For J1949+3 106, we 
used the Damour & Deruelle (henceforth DD) orbital model (IDamour fc Deruelldll986l ). with 
the Shapiro delay parameterized as in Freire & Wex (2010). The resulting best-fit parameters 
and other derived pulsar parameters are listed in Table [2j There are no significant trends 
in the timing residuals; this implies that at the current timing precision, the ephemerides 
presented in Table [2] describe the TOAs. Both pulsars have now been timed for four years. 
This has allowed very precise measurements of the spin and astrometric parameters, in 
particular precise determinations of the proper motions of these two new objects, to be 
discussed below. 



2.2. Polarimetry 

In order to measure the polarization characteristics of PSR J1949+3106, we observed it 
in coherent dedispersion mode at two different frequencies. Our 25 minute GBT observation 
used GUPP@ to sample a bandwidth of 200 MHz centered o n 820 MHz. We also observed 



the pulsar for 15 minutes at Arecibo using the ASP backend (IDemorestl 120071 ) with a band- 
width of 24 MHz centered on 1412 MHz. The data were analyzed in standard fashion with 
PSRCHIVE (Hotan et al. 2004), and the resulting calibrated full-Stokes pulse profiles are 
shown in Figure [T] and Figure HJ 



21 From Bureau International des Poids et Mesures (BIPM) 
22 https://wikio. nrao.edu/bin/view/CICADA/GUPPiUsersGuide 
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At 820 MHz, there is little to no circular polarization in the second component of the 
pulse profile, or linear polarization in both components. This is somewhat unusual since 
pulsars are among t he most polariz ed radio sources and individual pulsars can be up to 



100% polarized (e.g. lHan et al.ll2009f ). At 1412 MHz, there is some linear polarization in the 
first component, and the second seems unpolarized. The rotation measure is unconstrained 
in either observation. Comparing the two profiles, it is clear that the second component 
has a steeper spectrum than the first. This is confirmed by the observed profile evolution 
with frequency across the four WAPP bands used in Arecibo timing observations (Figure [3]). 
Assuming a 20% uncertainty, the measured period-averaged flux density is 0.39±0.08mJy for 
the 820 MHz GBT observation and 0.23 ±0.05 mJy for the 1412 MHz Arecibo obsservation, 
giving a spectral index of a = —0.97. 

We observed PSR J1955+2527 for 10 minutes using the same Arecibo ASP setup and 
data reduction method described above. The calibrated full-Stokes pulse profiles are shown 
in Figure HI There is little to no linear polarization and no circular polarization in the pulse 
profile, and the rotation measure is not constrained. The period-averaged flux density at 
1412 MHz is 0.28 ± 0.06 mJy. 



3. PSR J1955+2527 

PSR J1955+2527 is an isolated MSP with a period of 4.87 ms and dispersion measure 
of 209.97 pc cm -3 . Figure [3] shows a folded pulse profile from an Arecibo timing observation. 
Table [2] summarizes the timing solution parameters, and Figure [5] shows timing residuals vs. 
epoch for the TOAs after removing the best-fit timing solution. 

With its overall post-fit RMS timing residual of 11.4 fis, J1955+2527 does not have high 
enough timing precision to be included in the pulsar sample used by Pulsar Timing Array 
projects to attempt detection of gravitational waves^ 3 ! An RMS residual on the order of 
100 ns or lower is required for that purpose. In Section we investigate the contributions 
of various intrinsic and extrinsic effects to the timing residuals of both J1955+2527 and 
J1949+3106. 



2: 
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4. PSR J1949+3106 

PSR J1949+3106 is an MSP with a period of 13.14 ms and DM of 164.13 pc cm" 3 . 
It is in a 1.95-day binary system, and because the orbital plane of the binary is inclined 
by ~ 80° as viewed from Earth, we are able to measure the Shapiro delay in the system 
and derive estimates of the pulsar and companion masses, and the orbital inclination angle. 
Table [2] summarizes the timing solution parameters for J1949+3106, and Figure [6] shows 
timing residuals vs. orbital phase for the TOAs used to obtain the solution. Figures [Q [2j 
and[3]show folded pulse profiles for J1949+3106 at 820 MHz, 1412 MHz, and the four WAPP 
subbands, respectively. 



4.1. Keplerian parameters 

The timing solution of any pulsar binary includes the orbital period p,, the projected 
semi- major axis of the orbit x (typically expressed in light seconds), and parameters depend- 
ing on the orbital eccentricity e, the time of passage through periastron T , and the longitude 
of the ascending node cu. For PSR J1949+3106 they are presented in Table [2J From these, 
we can estimate the mass function of the system: 

„, . (m c sin i) 3 4n 2 x 3 c 3 / s 

f(m p ,m c ) = - ■ _ = ___, (1) 

(m p + m c ) <-r-r b 

where m p is the mass of the pulsar, m c is the companion mass, i is the orbital inclination 
angle with respect to the plane of the sky, c is the speed of light, and G is the gravitational 
constant. For PSR J1949+3106, we obtain / = 0.10938(5) M Q . Assuming that the pulsar 
has a mass close to the Chandrasekhar limit (1.4 M Q ) and an orbital inclination i = 90°, 
we obtain a minimum companion mass of 0.8 M . To disentangle the masses of the two 
components from the sine of the orbital inclination angle, we must be able to measure at 
least two post-Keplerian parameters ( 



4.2. Shapiro delay 

One important goal of pulsar surveys in general is to discover pulsars whose properties 
can constrain the equation of state of neutron-star matter. Current proposed equations of 
state differ in their allowed ranges for pulsar rotation periods and masses. Therefore, two 
ways of constraining them are to find very fast-spinning (P < 1 ms) or very massive neutron 
stars (m p > 1.8 — 2.0 M ). While pulsar rotation periods are obtained immediately upon 
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discovery from Fourier-based search algorithms, measuring pulsar masses is only possible for 
binary systems in which one or more post-Keplerian effects can be measured. 

When the companion is between the pulsar and the Earth, pulses traveling from the 
pulsar to the Earth pass through the gravitational well of the companion and a relativistic 
delay is intro duced in their arrival time as seen by the observer-the Shapiro delay (e.g. 



Shapiro! If 964J ). Most orbital models parametrize the Shapiro delay in terms of the "range" 



r = Gc 3 m c and "shape" s = sin %. Once we have fits for r and s, we can solve for m p based 
on the mass function of the system. 



Freire &: Wexl (120 101 ) rewrite the Shapiro delay in terms of two new post-Keplerian pa- 
rameters: h 3 , the amplitude of the third orbital harmonic, and the ratio of the amplitudes 
of subsequent harmonics. The correlation between the parameters of this "orthometric" 
model is much weaker than the correlation between r and s, leading to a better descrip- 
tion of the combinations of m c and sini allowed by the timing measurements. We use this 
variation on the DD binary orbital model for J1949+3106 (Table |2]). Figure [7] shows the 
constraints on the pulsar mass, companion mass, and orbital inclination that result from 
the measurement of <j and h%. It also shows alternative constraints derived using the normal 
r-s parameterization. The figure also shows 68.3% contours of the 2-D probability density 
function (PDF) derived from a x 2 map of the /13-/14 orthometric space, as described in sec- 
tion 5 of Freire & Wex (2010); this map provides very similar results to a x 2 map of the 
cosi-m c space. It is very clear that the 68.3% contours are better described by the s-h 3 
parameterization: there are points in the diagram (signaled in red) that are l-a consistent 
with r and s that provide bad (i.e., high x 2 ) fits to the timing data. 



4.3. Masses of the components of the J1949+3106 binary system 



The lateral panels of Figure [7] show the projection of the 2-D PDF into the m c , cosi 

M m and 



and m p axes. From these 1-D PDFs, we obtain m p = 1.47_q; 3 i M Q , m c 
i = 79.9^J'g°, where uncertainties correspond to ±l-er confidence levels. 



0.85; 



-0.14 

-0.11 



Intermedia te-mass puls ar binaries (IMBP) have C-0 or O-Ne-Mg WD companions with 
m c > 0.4 Mm jTauris 2011); for this reason we classify PSR J1949+3106 as a member of 
this class. Ferdman et al.l ( 2010 ) summarize possible IMBP evolution scenarios and argue 
that the IMBP J1802— 2124 has undergone common envelope (CE) evolution, leading to 
little mass accretion and a history similar to recycled pulsars with neutron star companions 



IMBP s with longer orbital periods tend to have larger eccentricities as well ( jTauris et al. 



20001 ) . J1949+3106 is in the orbital period regime where one can argue that CE evolution is 
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likely (Pb < 3 days). 



The observed eccentricities of IMBPs (e ~ 10 -4 — 10~ 5 ) are higher on average than 
those of low-mass pulsar binaries (LMBP), and the orbital eccentricity of PSR J1949+3106 
(4.3 x 10 -5 ) is well within this interval. Furthermore, the shorter accretion episode in IMBPs 
should not lead to the extreme quenching of the magnetic field seen in LMBPs. This also 
agrees with observa tion: IMBPs lie in a region of the P — P space distinct from LMBPs 
( ICamilo et al.ll20011r 4 l . with higher estimated magnetic fields (B oc \[P P). Again, the P, 
P and derived B for PSR J1949+3106 are consistent with those of previously determined 
IMBPs. 

Measuring the mass of PSR J1949+3106 precisely is important to test our understanding 
of stellar evolution. For IMBPs, the comparatively short accretion episode has another 
predictable consequence: t he total accreted m ass should be (much) smaller than the 0.1-0.5 
M Q expected for LMBPs ( jPfahl et al.l 120021 ). Therefore, the masses of pulsars in IMBPs 
should be very similar to their birth masses. This has recently been confirmed by th e mass 
measurement of the IMBP PSR J1802-2124 (1.24 ± 0.11 M Q ) ( TFerdman et~aDl201oh . This 
implies that for J1949+3106 we should also expect a mass close to the Chandrasekhar limit, 
and the current measurement of 1.471^31 M Q is consistent with this expectation. 

A precise measurement of three post-Keplerian parameters (<^, and Cj) will over- 
determine the mass equations for this system. This implies that this system will eventually 
provide a test of general relativity. 



4.4. Proper motion 



For pulsars with small period derivatives, the contribution of the Shklovskii effect 
( iShklovskiil Il970[ ) to P can be significant. The Shklovskii effect is due to the change in 
projected distance between the pulsar and the Solar System barycenter. Another effect 
contributing to the observed P is due to the difference in acceleration with respect to the 
Galactic center between the Sun and the pulsar. There is also a small contribution to the 
observed P due to the pulsar being accelerate d perpendicularly towards the Galactic plane 
( IDamour fc Taylor! Il99ll . iNice &: Taylor! 119951 ) . The contributions of these three effects to 
the measured P values for J1955+2527 and J1949+3106 are given in Table [3j The net value 
is subtracted from P listed in Table [2] for both pulsars before calculating the estimates for 



24 They note, however, that this is not consistent with classification based on m c : three of the systems 
have companion masses typical of LMBPs, 0.2 < m c < 0.3 M . 
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the spin-down luminosity E oc P/P 3 , surface magnetic field B oc and characteristic 

age r c = P/2P in the same table. 

One common feature of J1949+3106 and J1955+2527 is that the position angles of both 
pulsars' proper motions in Galactic coordinates are very closely aligned with the Galactic 
plane (G^ = 270° and 262°, respectively, Table [2]). A major reason for the proper motions 
being closely aligned with the Galactic plane is Galactic rotation. With the Sun and the 
pulsar moving around the Galactic center in different directions at ~ 220 km s" 1 , the relative 
velocity due to this motion is larger than the peculiar motion of the Sun and the motion of 
the MSP relative to the standard of rest at its position in the Galaxy. Another reason is 
a selection effect: the PALFA search region is in the Galactic plane, therefore MSPs with 
a significant component of motion away from the plane are less likely to be found. Given 
the current locations of the two pulsars (b = 2.55° and —1.58°, respectively), it is unlikely 
that they will ever move fa r from the Galacti c plane. In this respect they are ve ry similar to 
PSR B1855+09 (z ~ 50 pc. lKaspi et aJl994h and PSR J1903+0327 (z < 270 pc. lFreire et aL 



201 ll ) . which are also close to the Galactic plane. For the latter system, measurements of 
the radial velocity of the companion also suggest a relatively small velocity relative to the 
pulsar's LSR. 

Many ongoing s urveys are, for the fir s t time, finding a significant number of MSPs near 



the G alactic plane ( iCrawford et al.l 120121 . Keith et al.l l2010l iBates et al.l l201ll iKeith et al. 



201 ll ). Measuring the proper motions of these objects will be important to ascertain how 
many MSPs are tightly confined to the plane of the Galaxy. If there is a statistical excess of 
such systems then they could represent a separate, low-velocity MSP population. 



5. TO A Error Budget 

In this section we evaluate the relative contributions to the TOA precision that have 
so far been achieved for J1949+3106 and J1955+2527. First we test whether the residuals 
are consistent with white noise and then we evaluate the relative contributions to the white 
noise from three different effects. Given the sharp, unresolved main component in the pulse 
profile of J1949+3106, it is possible that broadband coherent dedispersion (e.g. with the new 
PUPPI backend, a GUPPI clone recently installed at Arecibo) might achieve significantly 
improved timing precision in the near future, providing a much more precise Shapiro delay, 
much more precise masses and a test of general relativity, as discussed above. 
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5.1. White Noise Tests 

A simple way to test w hether the po st-fit timing residuals are consistent with white noise 



is to count zero crossings ( jCordea 120111 ). This test is applicable to non- uniformly sampled 
data and is insensitive to discontinuities in the statistics of the white noise, for example 
jumps in variance due to using different instruments, as is the case with our residuals. For 
N samples of white noise, we expect on average (Z w ) = (N — 1) /2 zero crossings with a 
standard deviation of <Tz w = \JN — 1/2. A comparison between the observed and expected 
number of zero crossings for the residuals of J1949+3106 and J1955+2527 is shown in TableHJ 
including overall and per-observatory results. For 1949+3106, the number of zero crossings 
is within la of the expected white noise value if the data sets from each observatory are 
treated separately. For J1955+2527, the number of actual and expected zero crossings is 
within la for GBT TOAs and within 2a for all other observatories. 

Another white noise test involves fitting for a second frequency or period derivative. We 
expect that to be consistent with zero for white-noise-like residuals, and this is the result 
we obtain for J1949+3106. For J1955+2527, we detect a second frequency derivative of 
—2.2(7) s -3 , confirming the indication from the zero crossing tests that this pulsar's residuals 
may exhibit some red noise characteristics. 



5.2. Template Fitting Error 

Times of arrival are obtained by folding small portions of an observation (typically 
a few minutes) with an ephemeris that describes the timing solution of a pulsar, and then 
convolving the resulting folded profile with a pulse shape template. The pulse shape template 
used for extracting TOAs at a given observing frequency is the average of the folded profiles 
of many observations of the same pulsar at the same frequency. This is the method we use 
for obtaining TOAs for J1949+3106 and J1955+2527. 

Since TOA extraction is based on matched filtering, we can view the TOA and its 
uncertainty as a measurement of the pulse shape, which is affected both by the quality of 
the template and the amount of noise in the data. If the pulse shape does not change between 
observations, it ca n be described as an invar iant template added to noise. Assuming that 



the noise is white, ICordes fc Shannon! (120101 ) derive the minimum TOA error, a ts/N , from 



system parameters and an effective pulse width W e g. 



For a noiseless Gaussian pulse, W e ff = 0.6 FWHM. Using Equation A2 from lCordes fc Shannon 



( 120101 ) and smoothed versions of the TOA extraction templates for Arecibo observations of 



both pulsars with 1024 bins in the folded profile, we obtain W e g = 90 fis for J1949+3106 and 
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W cS = 109 fis for J1955+2527. For comparison, FWHM ~ 560 fis for J1955+2527, with un- 
resolved structure on the leading edge of the pulse deviating from a Gaussian. J1949+3106 
has a two-component main pulse with the FWHM of the brighter (and narrower) component 
of the main pulse ~ 180 /is. In both cases, W e s < 0.6 FWHM of the main pulse, because 
the TOA precision depends on the width of the narrowest substructure in the pulse that the 
template can match. In the case of J1949+3106, where the two components of the main pulse 
are well resolved, W e g ~ 0.6 FWHM of the narrower component. For J1955+2527, W e g is 
about a third of 0.6 FWHM. Even though the pulse profile of this pulsar is Gaussian-like, it 
is not completely featureless. The calculated value of W e g is likely affected by the presence 
of an unresolved bump on the leading edge of the main pulse and/or a slight bump at the 
very top of the pulse. 

For Arecibo L-band timing observations, T sys = 30 K and G = 10 K/Jy, giving S sys = 
3 Jy. The bandwidth per WAPP was 50 MHz, and each TOA was obtained by folding 
500 s of data. The peak pulse flux ~ 3.2 mJy for J1949 +3106 and ~ 1.8 mJy for 



J1955+2527. Using Equation Al from ICordes fc Shannon! ( 120101 ). we obtain cr ts/N = 4.6 /is 
for J1949+3106 and a is/N = 5.3 /is for J1955+2527. 



5.3. Pulse Jitter 



The folded pulse profile used in TOA extraction is obtained from averaging many pulses. 
For pulsars in general, each individual pulse is narrower than the average pulse profile, and 
the average pulse profile extends over the phase window where individual pulses are observed. 
Individual pulse phases may vary by an amount on the order of a pulse width from one pulse 
to the next. The folded pulse profile from which a TOA is extracted depends on both the 
shapes of individual pulses and the distribution of their phases. 

The intensity modulation index mi — eg/ (S) (the ratio of intensity rms and mean 
as a function of pulse phase) is used to chara cterize the amplitude modulation and the 



phase .fitter and is ty pic ally of order unity ( e. g. IHelfand fc Fowlerlll977l . iBartel et al.l Il980 



Weisberg et al. 1986 ). Cordes fc Shannon ( 201oh show that the TOA error due to pulse 
phase jitter, a tj , can be expressed in terms of mu the intrinsic pulse width W- m u arid a factor 



fj. Previous studies of profile stability (IHelfand et al.l Il975l iRathnasree fc Rankinl Il995h 
show results that are consistent with fj ~ 1/3 — 1/2 for most pulsars. Using Equation A6 
from ICordes fc Shannon! ( 120101 ). we calculate an upper limit on cr tj by assuming mj = 1, 
/j = 1/2, and W- in t equal to the FWHM of the 1.4 GHz template used to extract the Arecibo 
TOAs. With these parameters, we get cr tj = 1.2 fis for J1949+3106 and a t} = 0.52 /is for 
J1955+2527. 
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The average signal-to-noise ratio of an individual pulse is SNRi = SNR pro f / J N pu i scs , 
where SNR pro f is the signal-to-noise ratio of the folded TOA profile and iVp U i ses is the number 
of pulses averaged to produce the TOA profile. Generally, if SNRi > 1, the TOA uncertainty 
is dominated by pulse jitter as opposed to radiometer noise. Using our Arecibo TOA profiles, 
we obtain SNRi ~ 0.04 for J1949+3106 and SNR X ~ 0.01 for J1955+2527, consistent with 
the above result that a tj < <?t s/N and therefore the TOAs are radiometer noise dominated 
for both pulsars. 



5.4. Diffractive Scintillation 



The characteristics of scintles in any given portion of data that is p rocessed to produc e 
a TOA affect the folded pulse profile and introduce an error in the TOA (ICordes et al.lll990l ). 
We calculate this error for J1949+3106 and J1955+2527 from an estimate of the scattering 
broadening time r s at 1.4 GHz. 

An upper bound on the scattering time is simply W e «] therefore we have t s < 90 /zs for 
J1949+3106 and t s < 109 //s for J1955+252 7. For comparison, w e use the NE2001 model 
of the ionized gas distribution in the galaxy (ICordes &: Lazid 120031 ) for the coordinates and 
dispersion measures of the two pulsars to obtain r s at 1 GHz and scale it to 1.4 GHz assuming 
a Kolmogorov scattering spectrum with an index of —4.4. For an observing frequency of 
1.4 GHz, we obtain t s = 0.84 /is for J1949+3106 and 1.6 /is for J1955+2527, consistent 
with our upper li mits. As a further comparison, we use the empirical relation derived by 
Bhat et al.l (120041 ) . who assemble a set of pulsars with DM and r s measurements and fit 
a parabola with log(DM) as the independent variable. This relation gives t s = 8 /is for 
J1949+3106 and 127 fis for J1955+2527, though it is worth noting that there is scatter of 
up to two orders of magnitude in r s about the fit. 

The scintillation bandwidth A /diss an d time scale Ai^ss could not be constrained 
from the data s ince the two pulsar s are not bright enough to observe scintles. We find 
A fm ^ following lLambert &: Rickettl (119991 ) . We calculate At DISS using Equations 11 and 12 
from Cordes fc Rickettl ( 1998h and the di s tance D and transverse velocity Vr for each pulsar 
(Table [2]). Following ICordes fc Shannon! (120101 ). we calculate the number of scintles in time 
(N t ) and frequency (N{) for the portion of data used to produce each TOA. We obtain 
N t = 1.9 and N t = 51 for J1949+3106, and JV t = 2.2 and JV f = 101 for J1955+2527. Finally, 
we calculate the rms error in the scattering broa dening function due to scintillation, cr tsDISS , 
from Equation 23 in ICordes Shannon! (120101 ). Using the upper bounds on r s from the 
FWHM pulse width in the above calculations, we obtain cr tsDISS < 0.09 /is for J1949+3106 
and a tsDISS < 0.07 /is for J1955+2527. 
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5.5. Dispersion Measure Variations 

The dispersion measures of J1949+3106 and J1955+2527 and their uncertainties (Ta- 
ble^ were determined by TEMP02 fits to TOAs extracted from several frequency subbands 
per observation. While we have not observed DM variations in J1949+3106 and J1955+2527, 
we can place an upper limit on the TOA error due to unmodeled DM variations (whose effect 
could be partly or fully absorbed by the other parameters in the timing model) based on 
the DM uncertainties for both pulsars from the timing solutions in Table [21 The maximum 
contribution to the TOA uncertainty from unmodeled DM variations is 

_ 8.3/is <5DM Av 

.max rA \ / 

J center 

where Au is the bandwidth in MHz, /center is the center observing frequency in GHz, and 
5DM is twice the DM uncertainty. At the lowest frequency Arecibo observations used, 
1420 MHz, the bandwidth of 50 MHz implies DM uncertainties of 0.0006 and 0.003 pc cm" 3 
for J1949+3106 and J1955+2527, respectively. These values give cr tDMmax = 0.2 /is for 
J1949+3106 and 0.9 /is for J1955+2527. 



Backer et all (119931 ) fit observed gradients of pulsar DMs and find that the annual DM 



gradient is proportional to v 'DM. Based on their Figure 4a, we estimate an annual DM 
gradient of ~ 0.002 pc cm" 3 for J1949+3106 and ~ 0.003 pc cm" 3 for J1955+2527. For 
the ~ 4 years that we have been timing both pulsars, this empirical relation predicts a DM 
gradient contribution to the rms residual of 1.2 and 1.8 /is, respectively. Both predictions 
are within an order of magnitude of our estimates based on the DM uncertainties reported 
by TEMP02 above. 



6. Summary and Conclusions 

To the current tally of almost 2000 known pulsars we have added two MSPs found by the 
PALFA survey, J1949+3106 and J1955+2527, and presented their timing solutions. While 
J1955+2527 is isolated, in the J1949+3106 binary system we have been able to confidently 
measure the Shapiro delay and estimate the pulsar and companion masses. The pulsar's 
current median mass is 1.47 M Q , and the companion's median mass is 0.85 M Q . The uncer- 
tainties of these mass estimates are still 0.3 — 0.4 M and 0.1 — 0.2 M Q , respectively, but 
they will improve with continued timing. We are also on the verge of being able to detect 
the relativistic periastron advance in this system. 

We have outlined the steps towards breaking down the various contributions to the 
overall rms timing residual and applied them to the TOA sets of the two new pulsars. This 
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type of characterization is important with a view to figuring out what effects and in what 
cases are mostly responsible for the observed residuals, and what can we do to mitigate them. 
This has further implications for how we can take maximum advantage of the properties of 
new discoveries in projects like the International Pulsar Timing Array (IPTA), which need 
extremely precise timing measurements on pulsars that are very stable natural clocks. The 
North American Nanohertz Observatory for Gravitational Waves (NANOGra^y!^!) . an IPTA 
member, uses the Arecibo telescope for timing pulsars suitable for the IPTA and is especially 
interested in new MSP discoveries in the Arecibo sky. In the case of J1949+3106 and 
J1955+2527, the rms timing residuals are too large for including these pulsars in the IPTA 
sample. Table [1] lists overall rms timing residuals for both pulsars, residuals by observatory, 
the expected contribution to the residual from radiometer noise, and upper limits on the 
contributions from pulse jitter, diffractive scintillation, and unmodeled DM variations. We 
find that radiometer noise puts a hard limit on TOA precision in both cases. The overall 
timing residual of J1949+3106 is consistent with that limit, while the residual of J1955+2527 
is more than twice as large. The TOA residuals of J1949+3106 are consistent with white 
noise, while J1955+2527 exhibits some modest departures from white noise as evidenced in 
a zero-crossing test and in a fit for a second frequency derivative. 
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Fig. 1. — Polarimetric pulse profiles for PSR J19 49+3106 at 820 MHz. This is based on a 25 minute 
observation obtained at the GBT with GUPPI (jDemorest et al.ll2010i ). and is displayed with 256 
bins. In the bottom plot, the black trace corresponds to total intensity, while the red and blue lines 
correspond to linear and circular polarization, respectively. In the top plot, the position angle of 
linear polarization (PA) is plotted for bins in the the linear polarization profile with signal-to-noise 
ratio > 3. Bins outside the main pulse profile that exceed this threshold are likely due to statistical 
fluctuations or imperfectly cleaned RFI. 
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Fig. 2. — Polarimetric pulse profiles for PSR J1949+3106 at 1412 M Hz. This i s bas ed on a 



15 minute observation obtained at the Arecibo telescope with the ASP (jDemorestl 120071 ). and is 
displayed with 256 bins. In the bottom plot, the black trace corresponds to total intensity, while 
the red and blue lines correspond to linear and circular polarization, respectively. In the top plot, 
the position angle of linear polarization (PA) is plotted for bins in the the linear polarization profile 
with signal-to-noise ratio > 3. 
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Fig. 3. — Top: Pulse profiles for J1949+3106 from the four 50 MHz WAPP bands used in Arecibo 
timing observations. The two peaks have opposite relative strengths at L-band compared to 
820 MHz (Fig. [I]). Profile evolution is evident across the four WAPP bands: the trailing peak 
decreases in strength relative to the leading peak with increasing frequency. Bottom: Pulse profile 
for J1955+2527 from the 50 MHz WAPP band centered on 1420 MHz. Folded profiles from all 
Arecibo observations were aligned and averaged to produce these plots. For both pulsars, profiles 
use 1024 bins and have been normalized to the same height on the arbitrary vertical axis. The 
thickness of the line on the left middle of each plot corresponds to the width of a time bin. 
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Fig. 4. — Polarimetric pulse profiles for PSR J1955+2527 at 1412 MH z. This is based on a 
10 minute observation obtained at the Arecibo telescope with the ASP (lDemorestll20071 ) . and 
is displayed with 128 bins. In the bottom plot, the black trace corresponds to total intensity, 
while the red and blue lines correspond to linear and circular polarization, respectively. In 
the top plot, the position angle of linear polarization (PA) is plotted for bins in the the linear 
polarization profile with signal-to-noise ratio > 3. Bins outside the main pulse profile that 
exceed this threshold are likely due to statistical fluctuations or imperfectly cleaned RFI. 
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Fig. 5. — Timing residuals vs. epoch for PSR J1955+2527. Green TOAs are from GBT 
observations, blue TOAs are from Arecibo observations, cyan TOAs are from Jodrell obser- 
vations, and red TOAs are from Nangay coherently dedispersed data. 
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Fig. 6. — Timing residuals vs. orbital phase for PSR J1949+3106 before (top) and after 
(bottom) fitting for Shapiro delay parameters. The middle plot shows the part of the Shapiro 
delay that is not absorbed by the fits of Keplerian parameters of the system. Green TOAs 
are from GBT observations, blue TOAs are from Arecibo observations, cyan TOAs are from 
Jodrell Bank observations, and red TOAs are from Nangay coherently dedispersed data. 
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Fig. 7. — Constraints on the inclination angle and masses in the J1949+3106 binary system. The 
black contours include 68.3 and 95.4% of the total probability of a 2-D probability density function 
(PDF), calculated from a x 2 map of the h^-h^ orthometric space that used only the Shapiro delay 
to constrain the masses. The dashed purple lines indicate the constraints from /13 and the solid 
purple lines indicate the constraints from q. The solid black line indicates the upper limit derived 
for sjj.) The solid and dashed green lines show the constraints derived from the r-s parameterization 
of the Shapiro delay. The constraint based on s matches that derived from q, as expected. The 
r-s parameterization is sub-optimal: the points indicated by the red dots are 1-a consistent with 
the values of r and s, yet they provide very bad fits to the timing data. Nevertheless, both 
parameterizations provide very similar estimates of m c and cosi. left: cos i-m c plot. The gray 



region is excluded by the condition m p > 0. Right: 



m n -m, 



plot. The gray region is excluded by 



the condition sini < 1. Top and right marginal plots: 1-D PDFs for cosi, m p and m c , obtained by 
marginalization of the 2-D PDF. From the medians and ±1-<t percentiles of these 1-D PDFs, we 
derive m r 



1.47+°f M , m c 



0.85+°;}^ M and i 



79.9j£{j 
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Table 1. Contributions to TOA uncertainty (all values in /is). 



Parameter 


Observatory 


PSR J1949+3106 


PSR J1955+2527 


Timing residual 


Arecibo 


3.7 (1119) 


7.9 (199) 


(Number of TOAs) 


GBT 


12.8 (56) 


58.2 (97) 




Jodrell 


17.1 (85) 


36.2 (101) 




Nangay 


4.3 (28) 


11.8 (16) 


EFAC a 


Arecibo 


0.84 


0.96 




GBT 


1.09 


1.61 




Jodrell 


1.20 


1.22 




Nangay 


1.24 


5.70 


Overall rms timing residual 




3.96 


11.4 


Radiometer noise, Ats/N 


Arecibo 


4.6 


5.3 


Pulse jitter, Atj 


Arecibo 


< 1.2 


< 0.52 


Diffractive scintillation, Atsoiss 


Arecibo 


< 0.09 


< 0.07 


Unmodeled DM variations, At DM 








From DM uncertainty* 3 


Arecibo 


< 0.2 


< 0.9 


From DM gradient fit 


Arecibo 


< 1.2 


< 1.8 



a Weighting factor for TOA uncertainties chosen such that the ratio of reduced x 2 i n the 
TEMP02 fit to the number of degrees of freedom is ~ 1. 

b This paper, Table [2] 



Backer et al.l (119931 ) , Figure 4a 
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Table 2. Pulsar parameters. 



Parameter 



PSR J1949+3106 



PSR J1955+2527 



Right ascension (J2000) 
Declination (J2000) 
Proper motion in RA, fi a (mas/yr) 
Proper motion in DEC, /is (mas/yr) 
Spin frequency, / (s _1 ) 
Frequency derivative, / (s~ 2 ) 
Second frequency derivative, / (s~ 3 ) 
Epoch of timing solution (MJD) 
Dispersion measure, DM (pc cm -3 ) 



igh 49 m 29.63851 s (l) 
31° 06' 03.8289"(3) 

-2.94(6) 

-5.17(8) 
76.114023821963(3) 
-5.4407(3) x 10~ 16 

54500.000176077 
164.1264(5) 



19 h 55 m 59.39523 s (7) 
25° 27' 03.443"(2) 
-1.9(6) 
-2.4(8) 
205.22225531037(7) 
-3.84(4) x lO" 16 
-2.2(7) x 10~ 25 
54800.0 
209.971(3) 



Orbital period, Pb (days) 

Time of periastron passage, T (MJD) 

Projected semi-major axis, x (lt-s) a 

Longitude of periastron, oj 

Eccentricity 

Orthometric harmonic amplitude, h 3 
Orthometric harmonic ratio, q 



1.949535(2) 
54390.270(1) 
7.288650(1) 
207.5(2)° 
0.00004306(7) 
2.4(1) x 10~ 6 
0.84(2) 



Period, P (s) 

Period derivative, P (s s _1 ) 

Mass function (M ) 

Pulsar mass (M ) 

Companion mass (M Q ) 

Orbital inclination, % (deg) 

Rate of periastron advance, u (deg/yr) 



0.0131381833437039(5) 

9.3913(5) xlO^ 20 

0.10939(5) 
1 47+0.43 

-0.31 



0.85 
79.9 
< 0.02 



0.14 
-0.11 
1.6 
1.9 



0.004872765862979(2) 
9.12684(8) xl0~ 21 



Galactic longitude, I (deg) 

Galactic latitude, b (deg) 

Total proper motion, \x (mas/yr) 

Position angle of proper motion, © M (J2000) 

Position angle of proper motion, © M (Galactic) 

Surface magnetic field, B (G) b 

Spin-down luminosity, E (erg s _1 ) c 



66.86 
+2.55 
5.9(1) 
211(1) 
270(1) 
1.12 x 10 9 
1.63 x 10 33 



62.74 
-1.58 
3.1(8) 
218(13) 
276(13) 
2.14 x 10 s 
3.21 x 10 33 
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Table 2 — Continued 



Parameter 


PSR J1949+3106 


PSR J1955+2527 


Characteristic age, r c (Gyr) d 


2.2 


8.3 


Distance, D (kpc) e 


6.5 


7.5 


Transverse velocity, Vt (km s _1 ) 


~ 180 


~ 107 


Distance from Galactic plane, \z\, (kpc) f 


0.29 


0.21 


1400 MHz flux density, S1400 (mJy) 


0.23 ±0.05 


0.28 ±0.06 


1400 MHz radio luminosity, L 140 o (mJy kpc 2 ) g 


~ 9.7 


~ 15.8 


Number of points in timing fit 


1288 


413 


Weighted rms post-fit residual (fjs) 


3.96 


11.4 


Timing span (d) 


1715 


1773 



a x = a sin i/c, where a is the semi-major axis and i is the inclination angle. 
h B = 3.2 x 10 19 (PP) 1 ' 2 

C E = Att 2 IP/P 3 and assuming a 1.4 M Q neutron star with a 10 km radius and moment 
of inertia / = 10 45 g cm" 3 . 

d r c = P/2P 

c Esti mates based on DM, sk y position, and the NE2001 model of ionized gas in the 
Galaxy fjCordes fc Laziol 120031 ). 

f \z\ = D sin |b| 



g -^1400 — Si4Qo D 2 



Table 3. Contributions to measured P values. 



Effect 



PSR J1949±3106 PSR J1955±2527 



Shklovskii effect (s s _1 ) 


7(2) x 10- 


-21 


2.6(5) 


X 


10" 


-21 


Acceleration | to Galactic plane (s s _1 ) 


-6(1) x 10- 


-21 


-2.9(4) 


X 


io- 


-21 


Acceleration ± to Galactic plane (s s" 1 ) 


-7.8(7) x 10- 


-23 


-1.5(2) 


X 


io- 


-23 


Total P contribution (s s _1 ) 


0.274 x 10- 


-21 


-0.278 


X 


io- 


-21 
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Table 4. Number of TOAs, actual number of residual zero crossings, expected number of 
zero crossings for white-noise-like residuals, and expected standard deviation of the latter 
per observatory for J1949+3106 and J1955+2527. 



Observatory 


TOAs 


Z 


(z w ) 


°"z w 


J1949+3106 










Arecibo 


1119 


551 


559 


17 


GBT 


56 


31 


28 


4 


Jodrell Bank 


85 


42 


42 


5 


Nangay 


28 


16 


14 


3 


All 


1288 


613 


612 


17 


J 195 5 +2527 










Arecibo 


199 


109 


99 


7 


GBT 


97 


49 


48 


5 


Jodrell Bank 


101 


43 


50 


5 


Nancay 


16 


4 


8 


2 


All 


413 


206 


206 


10 
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